perm filename MODD.F4[JC1,MUS] blob sn#007329 filedate 1972-07-09 generic text, type T, neo UTF8
00100		DIMENSION  W1(512),FUNC(512),SU(50),ST(500)
00200		GO TO 60
00300	40	CALL HYDPOG(2)
00400	60	TYPE 50
00500	50	FORMAT(' TYPE A,CF,MF,I1,I2,STEPS'/)
00600		ACCEPT 100,A,FC,FM,XI1,XI2,STEPS
00700	100	FORMAT (6F)
00800		CALL TYPLOC(-300,-512)
00900		CALL DPYSET(1,SU,50)
01000		CALL DPYBRT(1)
01100		CALL AIVECT(0,0)
01200		CALL ALINE(-264,0,256,0)
01300		CALL ALINE(-256,-256,-256,256)
01400		CALL DPYOUT(1)
01600		XI=XI1
01700		XI1=(XI2-XI1)/STEPS
01800		N=STEPS+1
01900		DO 400 N1=1,N
02150		I=-1
02175		CALL DPYSET(2,ST,500)
02187		CALL DPYBRT(5)
02200		DO 200 J=1,512
02300		X=J-1
02400		Y=6.2832/512.
02500		FUNC(J)=A*(SIN((FC*X)*Y-XI*(COS((FM*X)*Y))+XI))
02600		IY2=FUNC(J)*256.
02700		IF(I.EQ.0)GO TO 150
02800		CALL AIVECT(-256,IY2)
02900		I=0
03000		GO TO 160
03100	150	CALL SVECT(1,IY2-IY)
03200	160	IY=IY2
03300	200	CONTINUE
03400		CALL DPYOUT(2)
03500		XI=XI+XI1
03600	400	CONTINUE
03650		TYPE 403
03675	403	FORMAT(' TYPE 0 TO CONTINUE'/)
03700		ACCEPT 300,M
03800	300	FORMAT (I)
03900		IF(I.EQ.0)GO TO 40
04000		END